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We present a detailed analysis of the contribution of small-angle Nambu-Goldstone boson 
(phonon) collisions to the shear viscosity, rj, in a superfluid atomic Fermi gas close to the unitarity 
limit. For low temperatures we show that the experimental values of the shear viscosity coefficient 
to entropy ratio, rj/s, can be reproduced assuming that rj is dominated by phonon collisions and 
considering an effective shear viscosity which takes into account the finite size of the experimental 
system. In particular, for very low temperatures, T < O.ITf, we find that phonons are ballistic 
and the shear viscosity is determined by the processes that take place at the interface between the 
superfluid and the normal phase. This result is independent of the detailed form of the phonon 
dispersion law and leads to two testable predictions: the shear viscosity should correlate with the 
size of the optical trap and it should decrease with decreasing temperature. For higher temperatures 
the detailed form of the phonon dispersion law becomes relevant and, within our model, we find 
that the experimental data for n/s can be reproduced assuming that phonons have an anomalous 
dispersion law. 

I. INTRODUCTION 

The study of the transport coefficients of a fluid opens a window on its microscopic dynamics, helping us to 
understand which are the underlying degrees of freedom and the most important interaction channels [l|-|3(. Of 
extraordinary interest is the study of the transport properties of fluids with an infinite two-body scattering length 
(corresponding to the so-called unitarity limit), because the expression of the transport coefficients in terms of the 
thermodynamic variables is believed to be universal [H, meaning that it is independent of the detailed form of the 
inter-particle potential. 

Experiments with trapped cold atomic gases are able to reach the region of infinite scattering length. When 
these experiments are performed with fcrmionic atoms, like 6 Li or 40 K, prepared in two different hypcrfinc states 
having a magnetic-field Feshbach resonance, the strength of the interaction between atoms with opposite spins can 
be varied by tuning the intensity of the magnetic field. At very low temperature, by varying the magnetic-field 
controlled interaction, fermionic pairing is observed to undergo the Bose-Einstein condensate (BEC) to Bardeen- 
Cooper-Schrieffer (BCS) transition, allowing us to explore the crossover region in a controlled way P-Q. Given 
the adjustability of the inter-particle interaction, ultracold fermionic systems can serve to emulate non-relativistic 
systems (which have the same spontaneously broken symmetries of the ultracold fermionic system) and to determine 
the universal expression of the transport coefficients at unitarity. On a more general ground, the results obtained can 
be used to emulate the properties of different systems which are believed to have similar properties as those of the 
Fermi gas at unitarity. One remarkable example is the fluid realized in rclativistic heavy-ion collisions, where almost 
ideal hydrodynamic flow has been measured [2|. 

Cold atomic systems at unitarity can also be used to explore the conjecture derived by string theory methods [ij], 
of the existence of a universal bound of the shear viscosity to entropy ratio rj/s > h/(4irkB)- The fluid produced in 
heavy- ion collisions seems to have extremely small shear viscosity close to the universal bound, while it is not yet 
clear which is the lower bound of r\js in fermionic cold atoms. As wc shall sec in detail, the present experimental 
results give values of rj/s above the universal bound, but the data seem to indicate that this quantity would further 
decrease with decreasing temperature. Recently, other different lower bounds on rj/s have been predicted based on 
the study of the effect of hydrodynamical fluctuations fiol . [Tl| . 

In the present paper we focus our attention to the superfluid phase of the ultracold Fermi liquid at unitarity. In 
superfluid Fermi systems the low energy degrees of freedom are the Nambu-Goldstone bosons (NGBs), which arise from 
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the quantum condensation of difermions by a process that spontaneously breaks a global continuous symmetry of the 
system. The property of superfluidity follows from the fact that the NGBs at low energy have a linear dispersion law. 
We will refer generically to these modes as supcrfiuid phonons, or phonons for simplicity. At very low temperatures 
they are the only relevant degrees of freedom, as the fermions are gapped and thermally suppressed. Thus, phonons 
give the leading thermal contribution to the thermodynamic and hydrodynamic properties of the ultracold superfluid, 
and the transport properties of the system strongly depend on the phonon self-interactions. 

Since the transport properties of the ultracold supcrfiuids are very sensitive to the phonon Lagrangian, it is possible 
to have hints on the phonon dispersion law and on the most important interaction channel by appropriate experimental 
measurements. In one of his pioneering works, Landau proposed a particular phonon spectrum (l2j, in order to explain 
the experimental values of the thermodynamic functions of 4 He. An analogous low energy spectrum can be inferred 
from the study of the transverse phonon relaxation processes, see e.g. [13[. Also the most relevant phonon self- 
interactions valid for 4 He can be determined in a similar way [Til [l5j . 

More recently it has been realized that the effective field theory (EFT) techniques can be used to determine 
the phonon dispersion law as well as the main phonon self- interactions for many superfluids [3]. The EFTs are 
very efficient when there is a hierarchy of widely separated energy or momentum scales (in our system, the phonon 
momentum versus the Fermi momentum), which allows one to obtain an expansion of the effective phonon Lagrangian 
in powers of energy or momentum, rather than in a coupling constant. Power counting and symmetry considerations 
fix the form of the EFT to the accuracy one desires, regardless of whether the underlying system is weakly or strongly 
coupled. At leading order in a momentum expansion, one can see that the EFT phonon Lagrangian is related to 
Landau's Hamiltonian by a Legendre transformation [l6|. Using the power of these EFT one can thus relate the 
physics of phonons in different supcrfiuid systems having the same global symmetries. 

In this manuscript we give an extended explanation of the results of Ref. , where we employed both the EFT 
techniques as well as the experimental available data, to pin the phonon Lagrangian down. In our procedure we 
assume a generic expression of the phonon effective Lagrangian, obtained by an expansion in powers of the phonon 
momentum, and study the contributions of various phonon processes to the shear viscosity coefficient of a superfluid 
ultracold Fermi liquid at unitarity. For the evaluation of 77, we consider the Beliaev 1 o 2 processes, using for this 
purpose the scattering rates obtained from the EFT Lagrangian of Ref. [lj| . We perform this calculation assuming 
an infinite volume and employing the results for the phonon spectrum obtained by various techniques [HQl]. 0ur 
results, and associated discussion, can be easily related to those obtained for 4 He in the temperature, T, regime 
where the viscosity is dominated by the phonons [Til [l5j . The resulting expression of the shear viscosity coefficient 
depends on one free parameter, which is related to the phonon dispersion law at high momentum and that has not 
been yet evaluated. Then, we try to fix this free parameter by comparing our results with those obtained by the Duke 
group (20I |22| , where measurements of the entropy of ultracold fermionic systems and of the shear viscosity coefficient 
(by measuring the damping of the breathing mode) have been obtained. Unfortunately reaching temperatures much 
smaller than the Fermi temperature, Tp, is extremely difficult and only few experimental points are available at the 
temperature where the contribution of phonons might be relevant. Moreover, the experimental setup is such that 
finite size effects must be included for an appropriate comparison with our results. Indeed, for T ~ 0.1 Tp, we find 
that for the particular experimental setup of Refs. [2(| phonons are ballistic, because their mean free path is of 
the same order of the size of the trapped cloud. Thus, the infinite volume limit we used in the evaluation of the 
shear viscosity coefficient and the whole hydrodynamical treatment of the phonons ceases to be valid. However, in the 
ballistic regime the collisions of phonons with the boundary still produce dissipation, and ultimately, the damping of 
the breathing mode. This dissipation can be described by an effective ballistic shear viscosity, employing an approach 
very si milar to the one used to describe the dissipation in the ballistic regime of the phonons of 4 He, or other systems 

Once the parameters of the effective Lagrangian of phonons are known, they can be used to predict the behavior 
of different transport properties at very low temperature, such as the thermal conductivity [26| , or the bulk viscosity 
coefficients [27| ■ In the present paper we content ourselves to consider the behavior of 77/ 1 s at very low temperature, 
marginally reached by today experiments, and to predict the behavior of this quantity with decreasing temperature. 
Indeed, according to our calculations, trapped phonons below ~ 0.1 Tp arc ballistic, and it follows that 77/s should be 
directly proportional to the temperature. We also find that rj should approximately scale with the size of the trap, a 
behavior that future experiments should be able to check as well. 

This paper is structured as follows. In Sec. |TT] we present a brief review of the phonon EFT at leading and next- 
to-leading order in the momentum expansion, with emphasis on the expression of the phonon dispersion law. We use 
the EFT to get the scattering rates necessary for the evaluation of the contribution to 77 of small-angle collisions in 
Sec. Ill Al In Scc. lIII Al we use a variational approach to solve the Boltzmann equation for three phonon collisions, and 
obtain their contribution to 77. The T dependence of such a contribution strongly depends on the form of the phonon 
dispersion law, and we obtain different values of 77 for some specific choices of the dispersion law in Sec. 1111 B[ In 
Sec.|lV]we discuss the experimental settings of the trapped Fermi gases where the measurements of 77 have been done. 
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In Sec. II V Al we compute the phonon mean free path and compare it with the size of the atomic cloud, to reach to the 
conclusion that finite size effects have to be taken into account. We write a phenomenological formula for 77 that takes 
into account those effects, and use it in Sec. IIVBI to give account of the experimental measures of r\j s. Appendix 
|A1 discusses why it is enough to keep the lower order corrections in the phonon dispersion law in the evaluation of 
77, while in Appendix [B] we present an alternative variational solution to the Boltzmann equation. Throughout the 
paper we use natural units K = ks = 1- 



II. PHONON COLLISIONS RELEVANT FOR THE VISCOSITY IN THE FERMI SYSTEM AT 

UNITARITY 

We call the superfluid phonon the NGB associated to the spontaneous symmetry breaking of a global U(l) symmetry 
group associated to particle number conservation. The phonon field can be viewed, by a microscopic approach, as 
the phase oscillation of the difermion condensate when the fermionic degrees of freedom have been integrated out. 
Actually, in this perspective, the low energy spectrum contains also the degree of freedom associated to the oscillation 
of the absolute value of the difermion condensate, see e.g. [28l - [30j . But considering only gaussian fluctuation, the 
radial degree of freedom can be integrated out and only results in a modification of the phonon dispersion law. In the 
following we assume that in the EFT Lagrangian the radial oscillations of the condensate have already been integrated 
out and the phonon is the only dynamical degree of freedom. It is worth to note that the fluctuations in the absolute 
value of the condensate are not taken into account in the analysis of [3J [HJ where the Lagrangian of the collective 
modes has been determined starting from a mean-field Hamiltonian. 

Effective field theory techniques can be used to derive the Lagrangian of the supcrfluid system at unitarity by a 
derivative expansion in the phonon field. At the lowest order, the effective Lagrangian has been obtained in Ref. |16| 
and is formally given at vanishing temperature by the functional 

£ LO = P(X) , (1) 

with 

X = (i -V(r)-d ti p-^L, (2) 

2m 

where P(p,o) and po are the pressure and chemical potential, respectively, of the non relativistic superfluid at T = 0, 
ip is the phonon field, and V is the trapping potential. For V = 0, the equation of state (EoS) of the unitary Fermi 
gas reads 

95/2 

rm 3/2 /4 /2 , (3) 



15tt 2 £ 3 /2 

where m is is the mass of the particles that condense and £ ~ 0.36 — 0.38 [33T - l3a | is the Bertsch number [36| . a 
universal constant that fixes the relation between chemical potential and Fermi energy, po = ^Ep- 

It is convenient to re-express the LO Lagrangian in a different way. After a derivative expansion, and rescaling of 
the phonon field to have a canonically normalized kinetic term, one obtains 

Clo = \ {{dt<t>) 2 - ^ h (V0) 2 ) - g {{d t <t>f - 3 5s &0(V0) 2 ) + Omf) + ■■■ . (4) 

where we have set V — for simplicity. The different self-coupling constants of Eq. (j4|) can be expressed in terms 
of derivatives of the pressure with respect to the chemical potential [37| . For comparison with the condensed matter 
literature on superfluidity, and for the purposes of computation, it turns out to be more convenient to express them 
in terms of the density p, the speed of sound at T = 0, c s , and derivatives of the speed of sound with respect to the 
particle density (see Appendix A of Ref. !37]) 




1 c 2 
v p h = c s = J — ^— , g= — — (l-2u) , 9s = - — V~ ' ( 5 ) 



1 - 2u 

where u = is the so-called Griineisen number. 

C s Op 

The next-to-leading order (NLO) Lagrangian is constructed by demanding non-relativistic general coordinate in- 
variance and conformal symmetry [ 1 €31 ] , and it reads 

£nlo = crm 1 / 2 !^! + ^(v 2 ^) 2 ^, (6) 
v X v TO 
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where c\ and ci are two dimensionlcss parameters, which are universal but cannot be determined by symmetry 
considerations. Both coefficients determine not only the corrections to the different phonon self-couplings in Eqs.©, 
but, and this is more important for the evaluation of the shear viscosity, also the corrections to the phonon dispersion 
law which at the NLO reads 

E k = c s k(l + 7 fc 2 ) , 7 = - (ci + |c 2 ) , (7) 



where hp = \/2mEF is the Fermi momentum. The sign of 7 has a dramatic effect on the possible phonon interaction 
channels, because, as we shall discuss in detail in the next section, the Beliaev process 1 <-> 2 between massless 
particles is only allowed for positive values of 7. The numerical values of c\ and C2 are still uncertain, but some 
progress has been made by a number of different techniques. Both coefficients can be estimated by the e-expansion 
technique, which gives Ci + 3/2c2 — —0.0209 [l8j|. However, these estimates include large errors, in particular for the 
determination of C2, because a finite value of ci only appears at the second order in the e expansion [l8j . A different 
approach is presented in (l9| where the phonon dispersion law is determined by a fit of Quantum Monte-Carlo (QMC) 
simulations. The reported result is c\ +3/2c2 ~ —0.028, but large uncertainties are present, arising from the numerical 
fit of the QMC data. A mean-field theory has been employed in Ref. [38j], giving c\ + 3/2c2 ~ —0.0205. Although 
with large numerical uncertainties all these methods give results consistent with a positive 7 coefficient. A different 
estimate of these parameters can be obtained by means of a mixed approach [39| . which combines the e-expansion 
technique and the results obtained considering the gaussian fluctuations on the top of the mean field solution f30| . 
According to [39[ , C2 is proportional to the square of the pair size and should therefore be positive (in disagreement 
with the estimate of [38[ ) and presumably comparable in size to C\. The results of [3!| lead to c\ +3/2c 2 ~ +0.01 and 
therefore to a negative value of 7, in disagreement with [ID, HH, [38|. However, considering gaussian fluctuations close 
to the unitarity might not be enough to determine the dispersion law of the collective modes, because the procedure 
seems to be not self-consistent [301 ■ For this reason we assume that 



0.18 

7 - ~W ' (8) 



1 F 

which is obtained considering c\ + 3/2c2 = —0.021 and £ = 0.37. We shall then briefly comment on how our results 
would change for negative values of 7. 

Given the expression of X in Eq.([2]), it is clear that in the NLO Lagrangian, not only the trapping potential appears, 
but also its space derivatives. In order to simplify the calculations we will use the local density approximation (LDA), 
which corresponds to consider the system as locally homogeneous, and therefore we will neglect the spatial derivatives 
of the trapping potential. Approximating the trapping potential by the harmonic potential 

V(x) = ±mu*x*, (9) 

the LDA turns out to be a good approximation as far as the local value of the chemical potential is much larger than 
u>i, which is a condition that is satisfied in the experimental settings of Ref. [2(|[22|, except near the edge of the cloud. 
In the following we will first consider the situation where there is no trapping potential, V = 0, and then comment 
on the effects in the phonon physics due to the presence of the trap. 

Terms in the phonon Lagrangian beyond the NLO will as well change the phonon dispersion law and the interaction 
vertices. We are not aware, to the best of our knowledge, of a determination of the phonon dispersion law at higher 
orders 1 . In general, the phonon dispersion law can be written as 

E k =c s k(l + i>(k)) , (10) 

where |^(fc)| <K 1 and it can be Taylor expanded in k 2 /k F , 

^ k ) =lk * + 6k A + o(jpj , (11) 

where S cx 1/kf, and we shall neglect terms of order higher or equal to k 6 /kp. We shall obtain a numerical estimate 
of 5 from a fit of the experimental values of the shear viscosity to entropy ratio. In doing this we shall follow a 
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procedure akin to the one used in 4 He, where the spectrum of the low energy excitations is determined from a fit of 
the experimental values of the mean free path of phonons, see e.g. (l3j . We comment on the appropriate temperature 
range of validity of the expansion in Eq. (fTTj) in the Appendix [Al where also an estimate of the k 6 /k F term is given. 

In principle, from a naive power counting, one would expect that NLO (or higher order) corrections to the interaction 
vertices should be taken into account as well. However, we shall show that in the computation of the shear viscosity 
coefficient these corrections are subleading with respect to the one determined by the dispersion law. The reason 
being that, as we shall discuss in detail in the next section, the dispersion law determines the power dependence of 
the shear viscosity coefficient on the temperature. On the other hand, the neglected terms in the interaction vertices 
do always give corrections of higher power in T/Tp. 

A. Three phonon and four phonon interactions 

For the computation of the shear viscosity coefficient one needs to evaluate the scattering rates of both the binary 
(hereafter 4-ph) collisions and the Bcliacv 1h2 (hereafter 3-ph) processes. The vertices associated to the correspond- 
ing interactions can be obtained by the momentum expansion of the Lagrangian in Eq.(TT]), and then the scattering 
amplitudes are readily obtained by standard techniques. 

Both the contributions of 3-ph and 4-ph processes to the shear viscosity depend on the dispersion law of the 
phonons. However, the 3-ph processes are more sensitive to the expression of ip( k ) because of kinematical restrictions. 
Considering the LO phonon dispersion law, E p = c s p, one finds that, for 3-ph processes, energy and momentum 
conservations only allow the collinear scattering, and the contribution of this process to the shear viscosity diverges. 
The reason being that the shear viscosity coefficient is proportional to the time scale for momentum transfer in 
the orthogonal direction to the flow, and if the collision is collinear it takes an infinite amount of time to transfer 
momentum in the orthogonal direction. 

Considering the generic corrections to the LO Lagrangian, the dispersion law has an additional term ip{ k )i see 
Eq. ([TUjl . which may favor or disfavor 1 <H- 2 processes depending on the sign of ip(k). If ip(k) < 0, 3-ph processes 
are kinematically forbidden. On the other hand, if ip(k) > (at least for some momenta k), 3-ph processes are 
kinematically allowed. In the following we shall assume that the NLO corrections give 7 > and therefore make i(>(k) 
positive for certain values of k. 

If the three phonons participating in the collision have momenta p, p' and k, then the collision angle between the 
phonons with momenta p and p' is given by 

cos v = 1 + A (pW) + k '^ k ') - (p + k ')i'(p + k ')) » (12) 

pi 

and since tp(k) <C 1, one can deduce that the deflection angle is small, and roughly proportional to 9 ~ ip 1 / 2 . After 

an analysis of the scattering processes one obtains that phonons with thermal momentum close to fc t h ~ 10T/c s give 

1/2 

the leading contribution to the viscosity, and one can estimate 8 t h ~ V'th ' wnere V'th = ipykth)- 

In the following analysis of the contribution of the 3-ph processes to the shear viscosity, we will be using a double 
expansion, both in the phonon momentum and in the collision angle. We will be considering the vertices of the 
scattering rates as arising from the lowest momentum expansion £lo j a s corrections from £nlo would be suppressed 
in the shear viscosity by powers of T/(c s kp) oc T/Tp. However, we will use the phonon dispersion law beyond LO, 
as the expression of ip(k) determines the collision angle, and approximate all different expressions in small angles. 
From the LO Lagrangian we have that the square of the scattering amplitude of the 3-ph processes is given by 

|X 3ph | 2 = 36 3 2 {-E p E p ,E k , + g a (E p p' • k' + E p ,p • k' + E k ,p ■ p')f , (13) 

where g s is defined in Eq. ([5|), E p , E p i and E^ are the energies of the phonons with momenta p, p' and k', respectively. 
After using the constraints of energy and momentum conservation and considering phonons with the dispersion law 
in Eq. ([TU1) . the LO scattering amplitude can be expressed as 

= 4cf(l + u) 2 2 
P 

This expression is formally equivalent to the scattering amplitude for the 3-ph processes in superfluid 4 He, and it has 
been used by Maris [l4[ and Benin [l5[ for the computation of the shear viscosity in that system. What changes are 
the explicit expressions of p and c s as well as the 0(ijj) corrections. The fact that the expressions of the 3-ph processes 
are the same in these two superfluid systems is not unexpected, as it was shown [l6| that the formal expression of 
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the effective Lagrangian at LO, Eq. ([T]), can be reiated by a Legendre transformation to the Hamiltonian proposed 
by Landau to derive the phonon self-interactions in superfluid 4 Hc. 

The expression of the scattering amplitude of the 4-ph process has been derived at the LO in momentum expansion 
in Ref. [40(. We refer the reader to that paper for the corresponding expression of |A^4 P h| 2 - 

The results obtained above still hold when including the trapping potential in the local density approximation. The 
phonon effective field theory, and the associated scattering amplitudes of the 3-ph and 4-ph processes, are modified 
by taking into account that the ^-dependence of the trapping potential can be included in a space dependent effective 
chemical potential, /i(x) = fXo — ^( x )- But this is not the only effect due to the presence of the trap. As we shall 
discuss in detail in Sec. IIV1 the restricted geometry dictated by the trapping potential leads to important boundary 
effects when the mean free path of phonons becomes comparable with the typical length scale of the trap. 



III. SHEAR VISCOSITY DUE TO SMALL-ANGLE COLLISIONS 



The shear viscosity rj measures the transport of momentum in the orthogonal direction to the hydrodynamic flow. It 
is generally dominated by large-angle collisions of the quasiparticle components of the fluid. There are cases, however, 
when small-angle collisions may dominate, as it may turn out to be more efficient to achieve a large-angle collision by 
the addition of many small-angle scatterings. In the low T regime of superfluid 4 He where the viscosity is dominated 
by phonons, the 4-ph large-angle collisions give the leading contribution to the shear viscosity only in a restricted 
range of temperature [4l[, while 3-ph small-angle processes dominate for any temperature below 0.7 K [T3 |. 

In this section we use a kinetic theory approach to compute the 3-ph contribution to the shear viscosity, assuming 
that the system is infinite and homogeneous. We use the scattering rate as arising from £lo evaluated in the 
previous section, but assume in the energy and momentum conservation the generic phonon dispersion in Eg. (1101) . 
with ip(k) > 0. We shall discuss in detail the range of validity of this method, that is why in the scattering rate the 
0(ijj(k)) corrections can be neglected. Eventually, we perform both an expansion in the phonon momentum and in 
the collision angle given in Eq. ([T^jl . finding that the temperature dependence of the shear viscosity strongly depends 
on the expression of the phonon dispersion law. 



A. Variational solutions of the Boltzmann equation 



When a shear stress is applied to a fluid, e.g. by contact with a sliding surface, it will perturb the fluid producing 
a gradient of the velocity field. In a superfluid, a sliding surface moving at sufficiently small velocity will only 
produce a gradient of the normal field velocity in the direction perpendicular to the surface. For small deviations 
from equilibrium, the energy-momentum tensor is given by 



STtj = -vVij = -v ( diV s + djVi - -<%V • V ) , (15) 



Tii = Z I T^s^FfM . (16) 



where V is the velocity of the normal fluid component. 

At the microscopic level, the phonon contribution to the energy-momentum tensor of the system can be computed 
by means of kinetic theory and it is expressed by 

^ 2 f d 3 p p t pj 
(2tt) 3 E p 

where f(p,x) is the phonon distribution function, which obeys the Boltzmann equation [l2| 

where we have assumed to be in the superfluid rest frame, and C[f] is the collision term. The processes that give 
the largest contribution to the collision term are 4-ph collisions and 3-ph splitting and joining processes. We shall 
assume that these processes can be evaluated separately (i.e. independently) and then compare the corresponding 
equilibration times. The contribution of the binary collisions to the shear viscosity has already been evaluated in 
Refs. 0, Therefore, here we shall focus on the 3-ph processes. For the 1 -s-> 2 processes the collision integral is 
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given by 



CW/] = ~ I \M^\\2*)H^\P - P> - K'){f{ P ){l + f{ P >)){\ + f{k>)) ~ f{p')f{k>){l + f{p))} 



+ I |^ 3 ph| 2 (27r) 4 ^ 4 )(P' - P K<) {f(p')(l + /(p))(l + J(fc')) - + /(?'))} , (18) 

Jp',fc' 

where we have denned the shorthand notation 

^ (19) 



(2tt) 3 2E p ' 

and the square of scattering amplitude A^3 P h was written in Eq. (|13|) . As expected, the collision term vanishes when 
evaluated using the phonon Bose-Einstein equilibrium distribution function, so that C[/ cq ] = 0. 

For the computation of the transport coefficients we consider small departures from equilibrium so that the dis- 
tribution function can be written as / = / cq + 5f, and then we linearize the transport equation in 5f. For the 
computations of the shear viscosity, one assumes that 

6f = -h(p) Pkl v J c * {1 + fc * ) , (20) 
where h(j>) is an unknown function, Vki is the tensor defined in Eq. (|15[) . and we have defined 

Pki = PkPi - T^hiP 2 • (21) 

Upon substituting the expression of the perturbed distribution function in Eq. ()16l) one can extract the value of 
the shear viscosity coefficient, which turns out to be 

^ = ^I^Wj^ + U)^)- (22) 
On the other hand, inserting Eq. (|20|) in the collision integral, we obtain at the linear order 

SC ee -±— Fl] [h(p)]V l3 , (23) 

and the linearized Boltzmann equation reads 

- /cq(1 + f ^ PijVij = -±- FlJ [h(p)]V l3 • (24) 



2pT r J J 2E P T 
With this last result one realizes that the shear viscosity can also be written as 



% P h = -^ j ' Pijh^F^hip)] , (25) 

or equivalcntly, 

% P h = ^= f (2ir)W\P -P>- A")|^ 3 ph| 2 /p cq (l + /;?)(! + f7) K) 2 - (26) 

01 Jp,p',k' 

where 

ee h(p) Pij - - h(k')k^ , (27) 

and we have used the property that the scattering amplitude is symmetric in the momenta of the three phonons. 

In order to have a lower bound on the shear viscosity coefficient it is possible to employ a variational method. In 
order to do so, one defines the inner product 



(28) 
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where \</>i) and |<fo) are two arbitrary functions. If we define the function \x) = Pij, then Eq. (|24|) can be schematically 
written as \x) = H\<f>), where $ is the solution of the linearized transport equation. The shear viscosity can be written 
as ?73ph = (xl^)j which is equivalent to Eg. (f22j) . or alternatively, 7^3 P h = as in Eg. (|26|) . Finally, from the 

Schwarz inequality (and considering that H is positive semidefinite) it follows that for an arbitrary function $, 

7731311 - WWW) ■ (29) 

which is saturated when $ is the exact solution of the linearized transport equation. 

We will compute the value of ?73 P h using this variational method, and assuming that the phonon dispersion law 
is of the form given in Eq. (jlOp. where ili{k) is treated as a small perturbation. In principle an upper bound of the 
shear viscosity can as well be obtained, e.g. by the method proposed in (43|. However, from the fact that ?73ph should 
diverge for tf>(k) = 0, it is possible to figure out the particular family of functions which maximizes Eq. (|29p . More in 
detail, we look for a solution that minimizes (A^-) 2 , and thus, considering the expression of H, maximizes the value 
of 773 P h- It is then easy to see that within the space of rational functions, the solution must be of the form 

h{p) = ^-(l + a^{p))+0^ 2 ) , (30) 

where a is the variational parameter. This form of the solution is such that (A^-) 2 oc tp 2 , while other forms of rational 
functions give a result of order one, and thus, bigger in the ^-expansion. We fix the value of a by minimizing 



A^ A^ = 1(42 + 24a + 6a 2 ) (pW) + fcV(fc') - (p' + k')^(p' + k')f , (31) 



which occurs when a = — 2. We have checked the stability of this variational solution by changing the form of 
variational functions used in the maximization procedure, see Appendix [B] We have also checked that our variational 
solution agrees with the variational solution of the phonon transport equation given in Ref. [l5j . where the shear 
viscosity due to joining and splitting processes for superfluid 4 He has been computed. 

With the form of h(p) reported above, we computed the shear viscosity coefficient at the leading order in tp, and 
we obtained 

?eV— — 

^ 15 7 cf M 

where 



^3ph=(-) + (3 2 ) 



M = / dp ' dk ' iP ' k ' iP ' + ^ bW) + " {P ' + k ' )W + ^ f ?'+ k ' {1 + fC p' )(1 + ] ' (33) 

and the star on top of the integral indicates that the integration is constrained by the energy and momentum 
conservation. As already remarked, the expression of the shear viscosity given above strongly depends on the form 
of the phonon dispersion law. We will explicitly check this fact in the following subsection by computing 77 3p h for 
different choices of ip(p). Notice that in Eq. (|33)) we did not need to include in the scattering amplitude the 0(ip) 
corrections. Actually, in Eq. (fl"4"]) we did not even calculate such corrections. We shall see at the end of the next 
section that neglecting these correction is equivalent to neglect temperature corrections of order 0(T/Tp) 2 in our 
final results. 



B. Values of the viscosity with different phonon dispersion laws 

Let us first consider a dispersion law such that tp(k) = -fk 2 (with 7 > 0), meaning that 3-ph processes are allowed for 
any value of the momentum of the interacting particles. By a power counting analysis one readily deduces the scaling 
rj oc 1/T 5 , which is the same scaling obtained with binary collisions of phonons [13, EH, and from the expressions 
derived above we obtain 

•'*^ 2 - lxl0 ~W§7^r (34) 
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Neglecting 0(ip) corrections, the entropy due to the phonons is given by 



Sph = W ' (35) 



and assuming the values £ = 0.37, and that |ci + 2>cij1\ ~ 0.021, we obtain 

^ ~ 4.0 x HT 9 5f . (36) 
s P h T s 

This should be compared with the shear viscosity coefficient obtained considering binary collisions of phonons. At 
the LO in the momentum expansion [l], |40| . |42| 

^^2.2.10- 7 H, (37) 

Sph -L 

which is about two orders of magnitude larger than 773 P h- Also notice that r?4 P h does not diverge for 7 = because no 
kinematical restriction applies to 4-ph processes. 

At this point it is important to remark that the above expressions lead to two unphysical results. Both ?73 P h and 
?74 p h diverge at vanishing temperature, and both arc below the experimental values (and below the universal limit) 
for sufficiently high temperatures. We shall see that these behaviors are due to an oversimplified description of the 
phonon fluid adopted so far, and they can actually be used to have information about the phonon fluid. We shall 
discuss the T — > limit in Sec. IIV1 and we shall see that the behavior of the shear viscosity in this limit is regulated 
by the geometrical extension of the trap. 

The decrease of rj^h/s with increasing temperature is due to the fact that a linear dispersion law was used in 
its computation (40l . |42| , while this behavior should be modified when considering higher order terms in the phonon 
dispersion law. Given that 7 cx 1/kp and treating 7T 2 /c 2 as a perturbative parameter, we expect that the first NLO 
correction to r/4 p h leads to 

(38) 



P h 



where B is some coefficient that depends linearly on 7, and thus vanishes for 7 = 0. The correction becomes more and 
more important when approaching T c and moreover with increasing temperature we expect that more terms should 
be included in the expansion of 7?4 P h in Eq. (|38|) . eventually preventing ?74ph/sph to decrease below the universal bound. 
Unfortunately, the evaluation of the coefficient B and of the other corrections is rather complicated. For this reason 
we shall not include in our discussion the contribution of 4-ph processes. 

The decrease of ?73 P h/s with increasing temperature is due to the fact we have assumed a dispersion law with 
"4>{k) = jk 2 . Since the typical momentum of phonons is of order 10 T/c s , by increasing the temperature, the neglected 
terms in the tp(k) expansion may become relevant. This motivates us to include terms of order k 4 /k F in the expansion 
of Eq. (|TT]) . that is, we take into account the next to next to leading order (NNLO) term in the momentum expansion 
of the dispersion law. While there are different computations in the literature of 7, the value of 6 has not been 
computed yet. We only know that 6 oc 1/kp. We assume that 6 < 0, as the opposite choice leads to even higher values 
of the viscosity at low T than those obtained in Eq. ([M)) . In the following section we shall determine the approximate 
value of S from a fit of the experimental data. 

With this choice of the sign of 5, the function ip is only positive in a restricted region of momenta. This has as a 
direct consequence that the three phonon processes are only possible when the phonon momenta are in a certain range 
of values. The constraints of energy and momentum conservation at the 3-ph vertex, implies that the integration 

ranges in the integral in Eq. (|33| are given by < k' < JA, where /3 = — and by p[ < p' < p' 2 , where 



p'1,2 = ( — i a/— 3fc' 2 + fg)/2. After defining the dimensionless parameter 



5/3 

t 
4(5 ' 



7 

Vw = -77 j (39) 



and 



2S T 2 

P = --~, (40) 

7 cj 
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we can express the shear viscosity coefficient as 



= VvT^+uy^ axI 0Y (41) 



where 

52 / j ;„./ i „.\\4/ e rfl/J , „,2 , _\\2 



/(/?) = /? 2 / dxdy(xy(x + y)) 4 (6 - 5/3(z 2 + y 2 + xy)) 2 f x+y (l + f x )(l + f y ) , (42) 

and x = p'c s /T and y = k'c s /T ', so that the integral is done over adimensional variables. This leads to the following 
expression of the shear viscosity to entropy ratio (we take £ = 0.37) 



V3 P h _ 0-09 1 Tp 
Sph ~ ^ax 10) T 4 



(43) 



The quantity (3 depends on T, and therefore the temperature dependence of r]3 P h/s is complicated. However, we 
find that for 0.003 < < 0.015, 10) « 2046.4, independent of T. In this region, we thus obtain that r/ 3ph cx 1/T, 
and then 

~ 4 " 2 2 10 ' ^l for 0.05 T F <T< 0.12 T F , (44) 

where the temperature range has been evaluated considering ip ms x ~ 0.2 — 0.4, see Sec. IIVBI 

For smaller values of the temperature we obtain that I {(3) cx /3 2 , leading to a viscosity 773 P h oc 1/T 5 , in agreement 
with the results obtained using a vanishing 8. This is as expected, because the leading contribution to 773 P h comes 
from phonons with momenta of the order 10T/c s , and at very low temperature, the correction introduced by the 5 
term of the dispersion law is of the order T 4 /(c s /cf) 4 ; and becomes negligible. We have not obtained an analytical 
expression of 10) for arbitrary values of (3, thus we will present only numerical results in Sec. HVI for particular values 

Of VAnax- 

From these results it is apparent that the 3-ph contribution to the shear viscosity coefficient is extremely sensitive 
to the dispersion law of phonons. At the LO (neglecting terms proportional to 7 and 5), the shear viscosity coefficient 
diverges (see the discussion after Eq. (fT3| ) . Including the NLO term proportional to 7 leads to r^ph cx 1/T 5 , and 
finally including the NNLO term proportional to S gives, in a certain temperature range, cx 1/T. These results 
have to be contrasted with the effect of the neglected terms in the vertices of interactions. These terms would lead to 
corrections to the shear viscosity coefficient proportional to (T/Tp) 2 (or higher powers) and are therefore negligible. 
In other words, the dispersion law determines the leading temperature dependence of ?73 P h, the vertex corrections lead 
to sub-leading temperature dependence. 

Given the strong dependence of the shear viscosity coefficient on the phonon dispersion law, one may question 
whether with increasing temperature higher order terms in the phonon Lagrangian may be relevant. In the Appendix 
[A] we argue that the neglected terms should be relevant only for temperatures T ~ (0.2 — Q.3)Tp; given that the 
transition to the normal phase takes place at T c ~ 0.16Tp 2 [3, [45|, those terms are irrelevant for a discussion of the 
shear viscosity coefficient in the supcrfluid phase. 



IV. SHEAR VISCOSITY IN TRAPPED FERMI GASES AND COMPARISON WITH THE 

EXPERIMENTAL VALUES 



The experimental measurement of the shear viscosity of ultracold fcrmionic atoms is performed confining a degen- 
erate mixture of an equal number of spin-l/2-up and spin-l/2-down in an optical trap. In the following we shall refer 
to the experimental setup realized in Refs. [2(| HH with TV ~ 2.0 x 10 5 6 Li atoms in a magnetic field of strength 
B = 840 G, close to the s-channel Feshbach resonance at B = 834.15 G [HI]. The unperturbed number density of the 
trapped fcrmions is given by 



Po(r) 



(2mE F f/ 2 



3/2 



3tt 2 



E 

i=i 



'2E F 



(45) 



2 This is the value of the critical temperature in an homogenous system 
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where m is the mass of the 6 Li atom and Ep is the Fermi energy of N free fermions in a harmonic oscillator potential. 
In Refs. [2(| H3, the transverse frequencies of the trap are u x ~ 2n x 1778 s -1 , uj y ~ 27r x 1617 s _1 and the net axial 
frequency is w z ~ 2-k x 76s _1 , corresponding to a Fermi temperature Tp ~ 2.4 10 _6 K. Given the asymmetry of the 
trap, it is convenient to redefine the radial coordinate as 

2 2 2 , 2 , 2 

w r =Wii + uj y y +lj z z (4oj 

where u> = (ajxiVyiVz) 1 ^ 3 . Using this definition of the radial coordinate, the system is spherically symmetric, indeed 
the trapping potential, defined in Eq. ©, is given by V(r) = ^uj 2 r 2 . The shear viscosity can be extracted at both 
high and low temperatures by studying either the anisotropic expansion of the atomic cloud or the breathing mode 
damping \It\ . 

The critical temperature of the trapped system is T c trap ~ (0.21 — 0.25)Tf [HI, H2|, and corresponds to the tempera- 
ture for which no region of the system is in the supcrfluid phase. Whether the system is locally in the superfluid phase 
depends on the local value of the Fermi temperature, and it is found that the system is supcrfluid for T < 0.16Tp(r). 
Here Tp(r) is the local value of the Fermi temperature in the LDA. The part of the fermionic cloud which is superfluid 
occupies the trap center, that extends up to a distance f(T) such that T = 0.167V(r). Since at very low temperature 
the matter distribution in the trap center is roughly constant, and the largest part of the fluid is in the trap center, we 
shall assume that Tp{r) = Tp for r < f. The trap center is surrounded by a corona of normal fluid, which corresponds 
to the range r > r, where T > 0.16Tp(r). For T > j£ ra P one has that r = 0; only fermionic modes are present and 
their collision processes determine the value of the shear viscosity coefficient. At high T a simple dimensional analysis, 
or a much careful study based on a Boltzmann equation, predicts the scaling rj oc T 3 / 2 [48J, in agreement with the 
experimental results [47J . For a recent determination of the shear viscosity using quantum Monte Carlo calculations 
see Refs. [H, H^. 

For temperatures below the superfluid phase transition, both the theoretical and experimental studies of the trans- 
port properties of the system are quite challenging. In this case the low energy spectrum of the system consists of 
two types of quasiparticles: the fermions and the phonons. The phonons are present in the center of the trap, where 
the supcrfluid phase is realized, while the fermionic modes are present in the trap center as well as in the corona. We 
shall assume that the leading contribution to the shear viscosity is due to the trap center, where the density is larger. 
In the trap center the system is supcrfluid and therefore the quasiparticle fermions are gapped, and their density is 
exponentially suppressed as e~ A / T , where A is the value of the fermionic gap. Nevertheless, their contribution to 



the shear viscosity may be sizable at any temperature. Indeed it was shown by Pethick et al. 51| that the shear 
viscosity coefficient of superfluid 3 He is not exponentially suppressed at low temperature, but tends to a constant value 
~ l/4r](T c ). The case of trapped ultracold fermionic system close to unitarity was studied in a recent publication (52[ 
and the contribution of fermionic modes to the shear viscosity was analyzed by computing the value of the effective 
carrier fcrmion number in the supcrfluid phase as a function of T. It has been concluded that the fermionic modes 
in the supcrfluid phase cannot be responsible for the experimental values of the shear viscosity below T c . Instead, it 
is suggested that the shear viscosity is dominated by the fermions in the corona of the trap. However, the obtained 
numerical values of rj/s are about a factor two larger than the experimentally measured values. 

Since phonons are gapless excitations, their number density is not exponentially suppressed at low T, but decreases 
as T 3 . Therefore, at sufficiently low temperature they should give a sizable contribution to the thermodynamic 
quantities as well as to the transport properties of ultracold fermionic systems. In the BEC limit it is indeed known 
that below T c , phonons give the dominant contribution to the entropy of the system as well to the shear viscosity 
coefficient fl2l |. On the other hand, in the BCS limit their contribution is negligible as compared to the fermionic 
contribution. At unitarity the contributions of fermionic and bosonic modes to the thermodynamic quantities might 
be comparable, as shown in [3fj| studying gaussian fluctuations on the top of a mean field solution. However, it is 
reasonable to expect that at sufficiently low temperature the contribution of fermions to the entropy of the system 
will be exponentially suppressed. Therefore, we shall assume that s = s p h, given by the expression in Eq. (|35D . 

Regarding the shear viscosity coefficient, one has to combine the fermionic and the bosonic contributions. One 
possible and challenging way of doing this, is to solve the Boltzmann equations which take into account both degrees 
of freedom. However, since the results of Ref. [52| show that the fermionic contribution to r/ is above the experimental 
values, there should be a mode with a lower shear viscosity coefficient. Indeed, r\ is directly proportional to the 
typical time for transport of momentum in the direction orthogonal to the flow, tj_, and therefore a mode with a 
shorter t± should exist. We assume that this mode is the phonon. Actually, the authors of Ref. [52|] claim that 
the interaction between phonons cannot lead to the experimentally measured value of the viscosity, because phonons 



do not couple to transverse probes. They also notice that the contribution to r\ found in Ref. [42j increases when 
the temperature is lowered, in disagreement with the experimental results. However, the arguments presented in 
Ref. [52| are not conclusive. As we shall show in the next section in the analysis of the low temperature regime, the 
inclusion of finite size effects in the phonon dynamics leads to a contribution to rj that decreases when the temperature 
decreases. Furthermore, there is no physical reason why a quasiparticle of the system that transports energy cannot 
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transport momentum in the direction orthogonal to the direction of the flow. In 4 He, indeed, it has been shown that 
at sufficiently low temperature phonons give a sizable contribution to the shear viscosity 12[ , and an analogous result 
should hold for any ultracold supcrfluid. 

When the temperature of the system approaches T^ rap from below, the contribution of fermions to the shear 
viscosity and to the entropy of the system should increase, for two reasons. In the first place, the radius of the trap 
center decreases, and a larger portion of the system is in the normal phase. Secondly, the value of the quasiparticle 
gap A should decrease with increasing temperature and therefore close to X'* rap the fcrmionic modes in the trap 
center should not be exponentially suppressed. Indeed, close to J"* ra P the thermodynamic description in terms of 
fermionic modes gives an excellent description of the experimentally measured properties of the system [2l], H|| . On 
the other hand, starting from vanishing temperature and increasing T, both effects mentioned above result in an 
increasing importance of the fermionic contributions, however it is not obvious at which temperature the fermionic 
contribution will dominate the bosonic one. At the present stage, only a comparison with experimental results can 
help to determine such a temperature. In the following we shall discuss the contribution of the bosonic modes to 
rj/s, and show that the very low temperature experimental points can be described by taking into account solely the 
contribution of phonons. 



A. Finite size effects 



As already mentioned in Sec. IIII Bl the contributions to the shear viscosity of 4-ph and 3-ph processes diverge at 
vanishing temperature. This unphysical behavior is due to the fact we have performed the various calculations in an 
infinite volume. In order to properly evaluate the contribution of phonons to the shear viscosity, it is important to 
consider the finite extension of the optical trap and the corresponding matter distribution. In the LDA the effect of 
the matter distribution can be taken into account computing the trap average of the shear viscosity evaluated in the 
previous section. At very low temperature, the matter distribution in the trap center is roughly constant, indeed we 
have assumed Tp(r) = Tp for r < r, and we expect that the trap average would lead to a tiny variation of our results. 

However, the finite extension of the trap does play an important role because the interaction of quasiparticles with 
the boundary can be neglected only if their mean free path is sufficiently small. More precisely, finite size effects have 
to be taken into account when the Knudsen number K n = defined as the ratio of the mean free path I of the 
quasiparticles over the typical size of the system L, is not small. 

The phonon mean free path Z p h can be obtained by considering the same 3-ph and 4-ph processes we have taken 
into account for the evaluation of the shear viscosity. The shear viscosity and the corresponding mean free path can 
be readily evaluated by the formula 

1 

f7buik = gPphC s i p h , (47) 

where p p h = 2 J 5< ^ is the phonon density of the normal fluid component. At low T, for the 3-ph processes one obtains 

that ?3 P h oc 1/T 5 , while for 4-ph collisions Z4 P h cx 1/T 9 . Thus, no matter which process dominates, for sufficiently low 
temperature the mean free path becomes larger than the typical length of the system. In the case we are interested in, 
we take as the typical size the smallest radius of the atomic cloud, R x , given in Eq. (|45)) . We present in FigQ]a plot 
of the value of l p h/R x as a function oiT/Tp for both 3-ph and 4-ph processes. The vertical dashed green line in this 
plot and in all the subsequent plots indicates the critical temperature T* rap = 0.25Tp reported in [22|. However, the 
critical temperature might be slightly lower, i.e. T^ rap ~ 0.212V, as obtained in [2l[. It has to be understood that the 
various plots of the phonon mean free paths are meaningful only for values of T < T* rap because phonons do not exist 
in the normal phase. For the same reason, any plot of the phonon shear viscosity coefficient is meaningful only in the 
superfluid phase. Regarding the 3-ph processes the mean free path strongly depends on the value of V'max- As we shall 
see in the next section, the experimental data seem to favor the values V'max ~ 0.2 — 0.3, and the lines named 3ph-a, 
3ph-b and 3ph-c correspond to the mean free paths obtained with V'max = 0.4, V'max = 0.3 and V'max = 0.2 respectively 
For V'max ^ 0.2 the mean free path of the 3-ph process is always larger than R x . For V'max > 0.3 the mean free path 
of the 3-ph process becomes of the order of the size of the atomic cloud only for T ~ 0.12 Tp, and rapidly increases 
for lower values of T. We also present the plot of hph/Rx, (named 4-ph) corresponding to the dashed black line in 
FigJT] Also in this case the mean free path becomes of the order of the trap radius R x for T ~ 0.12 Tp. Notice that 
the curves 3-ph-a and 4-ph almost overlap for T < 0.15Tp, however for larger temperatures 4-ph becomes extremely 
small, while in the 3-ph-a case it remains of the order of 0.1R X . Such a difference is due to the fact that the 3-ph shear 
viscosity has been evaluated including the term proportional to ip(k) in the dispersion law of phonons. In this case the 
phonon dispersion law has a maximum, which leads to a nonzero minimum of the corresponding mean free path for a 
certain value of T. The tp(k) term was not included in the evaluation of the 4-ph shear viscosity and for this reason, 
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the 4-ph curve of the mean free path does not have a minimum and should not be very reliable when getting close to 
ytrap _ a pj^jj we nave evaluated the mean free path for the 3-ph processes considering ?/> max 3> 1 (corresponding 
to 5 <C 7 2 , sec Eq. (|39|0 . meaning that the effect of the maximum in the phonon dispersion law manifests at very high 
momenta (and temperature). In the temperature range reported in Fig. [T]we have obtained a plot of the mean free 
path in concordance with the results that obtained for 5 — and very similar to the one of the 4-ph processes. 




FIG. 1: Plot of the Knudsen number, l p h/Rx, as a function of T/Tf for 3-phonon (3-ph) and 4-phonon (4-ph) processes. 
The black dashed line corresponds to the mean free path associated to the 4~ph process. The solid red line (named 3-ph-a) 
corresponds to the 3-ph process with -i/wx = 0.4, the dot-dashed blue line (named 3-ph-b) corresponds to the 3-ph process with 
VWx = 0.3 and the dotted red line (named 3-ph-c) corresponds to the 3-ph process with ^i m „ = 0.2. The horizontal black 
dashed line corresponds to l p h/Rx = 1. The hydrodynamic description is valid for l p h/Rx <S 1. The vertical dashed green line 
approximately corresponds to the transition temperature between the normal phase and the superfluid phase fH/ . 

When the mean free path of the phonons exceeds the typical size of the superfluid region, the transport properties 
of the phonons are mainly governed by their interactions with the boundary, rather than by the self-interactions. In 
the trapped atomic system the density is not constant, and the superfluid region at the trap center is surrounded by 
a corona of a dilute fluid of unpaired fcrmions. This layer effectively acts as a boundary for the superfluid phonons, 
which are confined to the superfluid center of the cloud, and when phonons reach this region they are absorbed or 
scattered back. If the process is elastic, then the phonon will only exert a normal pressure on the boundary, resulting 
in the propagation of low energy excitations in the normal phase corresponding to standard first sound. But, if the 
phonon is diffused, then it will as well exert a shear stress on the boundary inasmuch as the particles of a rarefied gas 
produce a shear stress when diffused by the walls of a container. 

In principle, both kind of processes, self-interactions and interactions with the boundary, take place for any value 
of the Knudsen number. However in the hydrodynamic regime, for K n <C 1, self-interactions arc the dominant 
microscopic processes leading to the largest contribution to the transport coefficients, while the interactions with the 
boundary are only relevant in a layer of length ^ p h (the Knudsen layer) close to the boundary. The interaction with 
the boundary can be effectively taken into account in the hydrodynamic equations by requiring that the fluid velocity 
vanishes at the boundary [23j-|25[, the so-called "no-slip" boundary condition. 

Finite size effects in the flow of rarefied gases and fluids have been studied thoroughly in the past (23l - l25| . Depending 
on the value of K n they can be taken care of either by imposing a "slip" boundary condition in the solution of the 
Navier-Stokes equations, or, at a microscopic level, by using the Boltzmann equation with appropriate boundary 
conditions. For values of K n < 0.1 (say) the Navier-Stokes equation are still valid, but the shear viscosity coefficient 
(and the transport coefficients in general) must be corrected including the finite size corrections, which for low Knudsen 
numbers can be expressed as 
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where ?7buik refers to the value of the viscosity in the hydrodynamic regime, and £ is a coefficient that characterizes the 
interaction of the quasiparticles with the boundaries for a particular geometry of the system. Eq. (|48[) gives account to 
a correction to the standard hydrodynamic equations proportional to K n . With increasing I, terms of higher order in 
K n must be added, but the computation of such terms is complicated even for an ideal gas [23l[24|. In any case, with 
increasing values of K n one is forced to abandon the Navicr-Stokcs equation and more moments of the distribution 
function have to be considered. For recent progress with Grad's moment method see e.g. [25j |. 

For values of K n >> 1, the system is in the ballistic (also named Knudscn) regime and no hydrodynamic description 
can be employed for the description of the phonon gas. In this case the interactions of phonons with the boundary 
are more frequent than the interactions among phonons, and the relevant relaxation time is the time between two 
collisions of a particle with the boundary, t&. As an admittedly rough model of the interaction between the phonon 
fluid and the fermions in the external corona, we assume that the boundary between the superfluid and the normal 
region is a sharp surface, i.e. a simple boundary [24 |. and, as in the Maxwell model (24l . |25I |. the scattering at the 
boundary is described by one single parameter, the accommodation coefficient x, meaning that particles impinging on 
the boundary are diffused with probability x and specular-reflected with probability 1 — X- This simplified description 
allows us to qualitatively capture the physical process happening at the interface, where a shear stress is exerted by 
the phonon gas on the fcrmionic fluid of the external corona. A more detailed and realistic modeling of the interaction 
of phonons with the normal fermions is certainly possible, however such a description would require the knowledge 
of the matter distribution, as well as the computation of the interaction of phonons with (almost) gapless fermionic 
modes. Moreover, for a proper description of the interface between the normal and superfluid component one should 
take into account the transfer of matter between the two phases. This effect can be taken into account by introducing 
one further phenomenological parameter, usually named the condensation coefficient (24j |. which is related to matter 
transfer between the two phases. We postpone a detailed analysis of the interface to future work and in the present 
paper we content ourselves with a phenomenological formula for the shear viscosity coefficient which has been used, 



derived and tested for the description of finite size effects in the flow of rarefied gases and fluids, see e.g. [53[, in the 
presence of a simple boundary. 

As in a rarefied gas, the phonons will exert on the boundary a shear stress which is proportional to the energy density 
of phonons and will damp the oscillation of the boundary. The same phenomenon takes place in 4 He at extremely 
low temperature, T < 0.5 K, where it has been observed that ballistic phonons can efficiently damp the movement 
of immersed objects. Experiments done with an oscillating sphere [5^ . a vibrating microsphere |55j . or a vibrating 
quartz tuning fork filj immersed in superfluid 4 He (see filj also for a comparison of the results obtained with various 
experimental apparatuses) show that the damping of the oscillations can be described, in perfect accommodation 
(x = 1)j by the introduction of a ballistic (effective) shear viscosity coefficient defined by 

r/b a u = \p v )xC s d, (49) 

where d is the typical size of the oscillating object. This expression leads to excellent agreement with the experimental 
data for large Knudsen numbers [4l|, predicting the correct dependence on the temperature and on the typical size 
of the oscillating object. 

Note that the expression above is equivalent to the one in Eq. (|4"T|) with the replacement of ^ p h with d and suggests 
that in the ballistic regime (and in the presence of a boundary) one can extend the standard definition of the transport 
coefficient by replacing the mean free path with the typical size of the system. Therefore, in our case, we define a 
ballistic (effective) shear viscosity as 

??ball = ^P P hCsX( c sn) = gP P hC s a, (50) 

where c s t\, = L is the typical size of the system, i.e. the extension of the superfluid region at the trap center (which 
is smaller than R x ). In the present paper we are not interested in a detailed description of the interface between the 
normal and the superfluid phase but rather to understand whether phonons can give a sizable contribution to the 
shear viscosity coefficient. For this reason, we define the phenomenological length a = xL, which we will use as a 
fitting parameter. 

In order to evaluate the shear viscosity in the intermediate region, K n ~ 1, one should employ the Boltzmann 
equation treating collisions of phonons with the boundary and collisions among phonons on an equal footing. Solving 
such a problem is complicated even for ideal gases [23l425| . However, since we know the behavior of the relaxation 
time for large values of K n as well as for K n -C 1, it is possible to employ the same reasoning at the basis of the 
Matthiessen's rule [56j | to obtain an expression of r] in the intermediate region, which interpolates between the values 
of the shear viscosity coefficient in the Knudsen and in the hydrodynamic regimes. For this purpose we define an 
effective relaxation time, r, incorporating the effects of inter-particle collisions and of the collisions with the boundary 

r-^r^+r^. (51) 
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This relation follows from the assumption that the total collision frequency is the sum of the frequencies of the 
two mentioned kinds of collisions, thus we are assuming that the two collision processes are not correlated. Since the 
shear viscosity coefficient is proportional to the collision time, we define the total effective shear viscosity as 



where ?7baii is defined in Eq. (|50[) and is the 3ph shear viscosity. In principle, the contribution of the 4ph collisions 
should be considered as well. However, the r^ph has been computed assuming ip(k) = 0. This approximation should 
not be of great impact in the evaluation of the shear viscosity at low temperature as shown in Rcfs. [4(| EH. However, 
it affects the higher temperature behavior of ?74ph, and for this reason it cannot be included consistently in our analysis. 
That is why we will neglect binary collisions in the following subsection. 



The experimental values of the shear to entropy ratio have been obtained by the Duke group in a number of different 
settings. In the present paper we refer to the data set reported in Ref. \at\ and we show that superfluid phonons can 
explain the low T experimental behavior of the shear viscosity coefficient. The caveat is that in our analysis we use 
several simplifying assumptions and two fitting parameters. Let us summarize our assumptions. We consider that the 
only contributions to the shear viscosity are due to processes taking place among phonons or between phonons and 
fermions at the interface between the superfluid and the normal phase. We neglect, instead, contributions to r] due 
to interactions between the fermions. In other words we assume that any oscillation of the corona is damped by the 
friction with the trap center, and this process is effectively described by the ballistic shear viscosity in Eq. (jST))) . With 
increasing temperature, the contribution to ?y of the interaction processes among fermions should become more and 
more relevant due to the fact that the spatial extension and the density of the external corona increases and because 
more Bogolyubov modes are available. Therefore our results should be reliable only up to some temperature T < T c . 

Regarding the phonon contribution, we also neglect corrections to the shear viscosity coefficient coming from the 
higher energy part of the phonon spectrum. This amounts to neglect terms of order (k/kp) 7 in the phonon dispersion 
law. In the Appendix [A] we show that these terms should be negligible for temperatures up to (0.2 — 0.3)TV. 

While our formalism allows us to compute the shear viscosity in any point of the trap, performing the trap average 
is not a simple task. Thus, we assume that the largest contribution comes from the trap center, where the density is 
larger and almost constant, and where the local Fermi temperature, Tp(r), is roughly constant and we approximate 
it with the global value Tp. 

In order to explain the low T data we will use Eq. (|52[) . which depends on two unknown parameters, a and 5 (or 
equivalently i/) max ). We remind that the value of S could be in principle determined if the phonon dispersion law is 
computed to next to next to leading order, while a requires to know both the size of the superfluid core in the trap 
and the probability of the phonons to be diffused at the edges of the core. 

In Fig.[2]the experimental values of rj/s obtained in [47} are shown, together with three fits obtained by considering 
the contribution of the interaction of phonons with the boundary, described by Eq.(|50p. assuming a = 0.3 R x , and 
the 3-ph processes with ip max = 0.2 (dashed black line), if> max = 0.3 (dotted red line) and t/Wx = 0.4 (solid blue 
line). The experimental data seem to favor the value f/'max = 0.3. However, more experimental data are needed in 
order to figure out the correct value of V'max- Note also that the experimental values of rj/ s decreases with decreasing 
temperature and this behavior is well reproduced for i/) max < 0.3. Although all the experimental values arc above the 
value 1/4tt it seems that further decreasing the temperature values below the universal bound might be reached. 

In Fig. [3] we present three fits of the experimental data considering the contribution of the 3-ph processes (with 
"0max = 0.3) and the interaction of phonons with the boundary, effectively describe by Eq. (|50")) with a = 0.2 R x 
(dashed black line), a = 0.3 R x (dotted red line) and a = 0.4 R x (solid blue line). The experimental data seem to 
favor a = 0.3 R x , but more experimental data are needed to determine the precise value of this parameter. 

In Fig[?]we report several plots of all the various contributions separately as well as the effective shear viscosity. The 
interaction of the phonons with the boundary (dashed black line) is obtained from Eq. (|50[) assuming a = 0.3 R x ; the 
contribution of the 3-ph process (dotted red line) is obtained from Eq. (|43)) considering ip max = 0.3; the contribution 
of the 4-ph process (dot dashed green line) is obtained from Eq. (|3T|) . The effective shear viscosity (solid blue line) is 
obtained combining ^baii and 7?3 P h in Eq. (|52p with the above reported values of a and f/'max- With this choice of the 
parameters the experimental values are appropriately reproduced. 

We refrain from showing results for the effective shear viscosity including 4ph processes. Actually, for sufficiently 
low temperature the expression of 7?4 P h reported in Eq. Q37|) is reliable and it can be included in the computation of 
the effective shear viscosity; but in this case the contribution of 4ph processes is however negligible, because of the 




(52) 



B. Comparison with the experimental data 



16 




FIG. 2: Plot ofn/s in units ofh/ks as a function of T /Tf considering the contribution of the interaction with the boundary, 
assuming a = 0.3 R x , and the 3-ph processes with three different values of ^ maj . The dashed black line is obtained with 
f max = 0.2, the dotted red line is obtained with ip max = 0.3 and the solid blue line is obtained with i/Wx = 0.4. The experimental 
values and the corresponding error bars were taken from fTftl . The vertical dashed green line a ppr oximately corresponds to the 
transition temperature between the normal phase and the superfluid phase in a trap reported pZzll . The actual value of T* rap 
might be slightly lower; the value reported in is T* rap ~ 0.21Tf- 



dominant role of the ballistic viscosity. For higher temperature Eq. ([57)) is not reliable and NLO contributions of the 
phonon dispersion law should be included. 

Our plots show that at very low temperatures, below 0.1 Tp, the ballistic shear viscosity dominates, meaning that 
the only relevant dissipative mechanism is the one that takes place at the interface between the normal and the 
superfluid phase. This result is independent of the detailed expression of the phonon dispersion law, in particular 
this result remains correct also if 3-ph processes are not allowed, i.e. considering 7 < 0, and only binary collisions are 
included in the analysis. It only depends on the fact that at sufficiently low temperature phonons become ballistic, 
see Fig. [TJ 



V. SUMMARY 



Let us give a very brief summary of our results. Below the critical temperature the trapped Fermi gas is in a 
superfluid phase, characterized by the spontaneous breaking of a U(l) symmetry. In this regime there are two types 
of quasiparticles: the fermions and the phonon (Nambu-Goldstone) modes. At sufficiently low T, only the phonons 
give relevant contributions to the thermodynamics and transport coefficients of the system, as the fermions are gapped 
and thermally suppressed. At unitarity, EFT techniques and numerical methods can be used to write down the phonon 
Lagrangian and to assess the main phonon self-interactions, but this information it still incomplete, because not all 
the terms in the Lagrangian beyond the NLO have been determined, so far. Further and complementary information 
can be deduced by the analysis of the experimental data. This mixed approach, already used for different superfluids 
such as 4 He, allows us to constrain some parameters of the effective low energy Lagrangian and to understand the 
most important physical processes. 

We have presented a detailed computation of the contribution of small-angle collisions to the shear viscosity. The 
T dependence of such a contribution strongly depends on the form of the phonon dispersion law, as it had previously 
been discovered in 4 He [l4| • In this regard, one notices that the knowledge of the value of r\ in the T regime where it 
is dominated by phonons would certainly provide detailed information on their dispersion law. 

We have then discussed the relevance of the restricted geometry for trapped Fermi gases in the evaluation of the 
shear viscosity. Below T* Tap the superfluid phase is realized in the core of the gas cloud, while being surrounded by 
an outer layer of fermions in the normal phase. At T ~ 0.12TV, the phonon mean free path becomes of the same 
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T/T p 

FIG. 3: Plot of rj/s in units ofhjks as a function o/T/Tf considering the contribution of the interaction with the boundary and 
the 3-ph processes with ^ max = 0.3, and a = 0.2 Rx (dashed black line), a = 0.3 R x (dotted red line) and a = 0.4 R x (solid blue 
line). The experimental values were taken from XI\] . The vertical dashed green line approximately corresponds to the transition 
temperature between the normal phase and the superfluid phase EH/ . 

order, or even larger, than the size of the cloud, and phonons are then in the ballistic regime. Then, phonons collide 
more often with the boundary of the superfluid core, that is, with the outer layer of normal fermions, than among 
themselves. Such collisions can provide a shear stress on the boundary, which then results in the damping of the 
breathing modes. This dissipative mechanism can still be effectively described by a ballistic shear viscosity. The same 
phenomena occur in He, where it has been observed that ballistic phonons can damp the movement of immersed 
objects in the superfluid. 

In order to describe the measured values of r\ for temperatures in between the ballistic regime and the hydrody- 
namic regime, we have used a phenomcnological formula for 77, given in Eq. (|52p . which contains two (yet) unknown 
parameters, that we fit to the data. One parameter describes the phonon dispersion law at ncxt-to-next-to-lcading 
order, which in principle could be computed from the microscopic physics, while the other one is related to the size of 
the superfluid core in the trap and the type of phonon-fcrmion scattering processes that take place at the boundary. 

Although our model relies on a number of simplifying assumption, detailed in the previous section, and, as a 
consequence, gives a rough description of the experimental system, it nevertheless allows us to make a couple of 
testable semi-quantitative predictions. First, we notice that if experiments are conducted reducing/increasing the 
size of the trap, but keeping Ef constant, then rj/s at low temperatures should decrease/increase. In other words, 
the value of rj/s should correlate with the size of the gas cloud. Second, we predict that rj/s should decrease with 
decreasing temperature. If we naively extrapolate our results to lower temperature, we predict that there should be a 
violation of the string theory proposed bound of rj/s. Note, however, that this violation happens because phonons are 
ballistic, while the string theory bound concerns the hydrodynamic regime. Both predictions are independent of the 
detailed form of the phonon dispersion law, in particular, they are independent of the sign of the 7 term in Eq ([7|). 
However, if 7 > 0, meaning that in a certain range of momenta the Beliaev processes are allowed, it is possible to 
deduce further information on the phonon dispersion law from the experimental data. In particular, we find that the 
correction to the dispersion law should be given by the expression reported in Eqs. (fTO]) an (fTTj) with 5 ~ 0.03/fc^. 
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FIG. 4: Plot ofri/s in units of h/ks as a function o}T/Tf- We report all the various contributions separately as well as 
the effective shear viscosity obtained with Ea. \§2\ . The interaction of the phonons with the boundary (dashed black line) is 
obtained from Eg. (|50|l assuming a — 0.3 R x ; the contribution of the 3-ph process (dotted red line) is obtained from Eg. (|43[) 
considering V'max = 0.3; the contribution of the J^-ph process (dot dashed green line) is obtained from Ea. $i7\ . The experimental 
values and error bars were taken from Tvh l. The solid blue line is obtained with Eq. (|52l) combining the ballistic term and 3-ph 
processes. The vertical dashed green line approximately corresponds to the transition temperature between the normal phase and 
the superfluid phase fil/ . The actual critical temperature may be lower, T* rap ~ 0.212> as reported in 



Programme from Ministerio de Economi'a y Competitividad and from FP7-PEOPLE-2011-CIG under Contract No. 
PCIG09-GA- 2011-291679. 



Appendix A: Phonon dispersion law and possible "roton" excitations 

From the fit of the experimental values of the shear viscosity performed in Sec. IIVBI we have determined the 
coefficient S of the phonon dispersion law in Eq. (jlip . In order to emphasize the region of validity of the expansion of 
the phonon dispersion law obtained, we rewrite it as an expansion in x = k/kF, 

E/c s k F = x ( 1 + yx 2 - -^xA , (Al) 

where we have used Eq. ([M| . and we have defined = j/kp. The corresponding plot, taking ?/> max = 0.1 (solid black 
line), -0max = 0.2 (dashed blue line), ipm ax = 0.3 (dotted red line) is reported in Fig. [5] for positive values of E. 

This plot allows us to figure out whether the neglected terms in the phonon dispersion law are relevant for the 
present analysis of the shear viscosity coefficient. The terms not included in the momentum expansion are of the order 
of x 7 — (fc/fcp) 7 and one would naively expect that they should give a large contribution to the dispersion law, and 
to the integral in Eq. (J33J) , as soon as k ~ kp. However, the expansion in Eq. (|Al[) suggests that the actual expansion 
parameter is not 2, but i/jx and therefore the naive power counting in powers of k/kp is incorrect. Assuming that 
the expansion is in powers of \f^x and considering that 7 ~ 0.18, sec Eq.([7|), it follows that Eq. (|A1[) gives a good 
approximation as far as k < 2kp. This result is corroborated by the observation that the neglected terms should 
lead to an increase of the phonon dispersion law at k ~ (2 — 3)/cf in order to avoid that the dispersion law becomes 
negative, but should not be relevant at smaller values of k (otherwise it would lead to the wrong behavior of the shear 
viscosity at low temperature). Therefore a term, presumably of order ^(k/kp) 7 should exist, with positive sign, and 
lead to a minimum of the dispersion law; such a minimum should correspond to rotons. Regarding the temperature 
range at which such a term would be relevant, we can estimate it considering that the largest contribution to the shear 
viscosity coefficient comes from momenta k ~ 10 T/c s . Thus the present analysis is consistent up to temperatures 
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FIG. 5: Plot of the dispersion law of phonons given in Eq. (|A1|) . with tp n 
line), i/Wx — 0.3 (dotted red line). 



0.1 (solid black line), ip n 



0.2 (dashed blue 



T ~ (0.2 — 0.3) Tp, that is for temperatures of the order or larger than the transition temperature. Thus, the neglected 
part of the spectrum should be important only close to T c . 

We note that the presence of a dip in the phonon spectrum, both in the BCS and BEC regimes, has been excluded 
by the study of gaussian quantum fluctuations on the top of a mean field solution (30j . However, the procedure 
employed in this calculation is only reliable in the BCS and in the BEC limits and it is not clear how a self-consistent 
calculation might be done close to unitarity (see Section IX of [3(|). Moreover, higher order fluctuations of the radial 
oscillation of the condensate (not included in that analysis) might be relevant close to unitarity, possibly affecting the 
phonon dispersion law. 

Appendix B: Variational treatment with orthogonal polynomials 

In this Appendix we show that with a different choice of family of variational functions we reach to the same value 
of the shear viscosity as that one found in Sec. IIII Al Following a treatment very similar to the one carried out in 
Ref. 143 we select the trial functions as 



h(p)=p n Y,b s B s (p) , 



(Bl) 



where n is a parameter that will be determined by a variational procedure, and B s (p) are orthogonal polynomials of 
order s and are defined such that the coefficient of the highest power p s is 1, and the orthogonality condition 



d? P /eqCl+Z eq). 

(2tt) 3 2E p 



-Pijpijp 



l B r (p)B s (p)=A^6 rs 



(B2) 



is satisfied. 

Using this form of the solution one can check that Eq. (|22|) gives 



2c t, An) 



(B3) 
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where 



4 n) = ^5^r(6 + n)C(5 + n), (B4) 



where T(z) and £(z) stand for the Gamma and Ricmann zeta functions, respectively. 
On the other hand, the expression of the shear viscosity given by Eq. (|2"5|) reads 



JV=oo 



r?3ph = J2 h ^ M ^ > ( B5 ) 



s,t=0 



where 



and 



M st = -L / (2^) 4 <5( 4 )(P - P> - K')\M\ 2 r p \l + fP)(l + f^Atj 

01 Jp,p',k' 



b = ^A n \M 



so that 

4cf 



25T 2 



(B6) 



A\ 3 = B t (p)p n Pij - B t { P ')p' n p' l3 - Bt {k')k ln k' l3 , (B7) 
Requiring that the two forms of the shear viscosity, Eqs. (|B3j) and ()B5j) . to be equal implies 



(A^f(M- l ) 00 . (B9) 



In practical terms one performs the study by limiting the number of orthogonal polynomials included in the analysis. 
One can prove that 

^>^{A^f{M- l ) QQ , (BIO) 

for a particular value of n and the number of orthogonal polynomials considered in the study. 

For definiteness, we consider ip(p) = jp 2 - We then find that n = —1 leads to the maximum value of 773 P h- Further, 
after using only up to the polynomial of order s = 2we get the convergence to an expression which matches that of 
Eq. ([MD. 
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